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Possibility of a ferromagnetic semiconductor single wall carbon nanotube (SWCNT) , where ferro- 
magnetism is due to coupling between doped magnetic impurity on a zigzag SWCNT and electrons 
spin, is investigate. We found, in the weak impurity-spin couplings, at low impurity concentrations 
the spin up electrons density of states remain semiconductor while the spin down electrons density 
of states shows a metallic behavior. By increasing impurity concentrations the semiconducting gap 
of spin up electrons in the density of states is closed, hence a semiconductor to metallic phase tran- 
sition is take place. In contrast, for the case of strong coupling, spin up electrons density of states 
remain semiconductor and spin down electron has metallic behavior. Also by increasing impurity 
spin magnitude, the semiconducting gap of spin up electrons is increased. 



I. INTRODUCTION 



II. MODEL AND FORMALISM 
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Single wall carbon nanotubcs extensively attracted 
from both theoretical and experimental point of view, 
due to it's technological applications such as nano elec- 
tronics devises^, nano ring (SWCNT ring)£, molecule- 
nanotube hybrids^, quantum dots^. Effect of substitu- 
tion non magnetic impurity on a SWCNT are investi- 
gated by many people o 6l7 i- . We found that by dop- 
ing boron and nitrogen the semiconductor gap, E g , of 
a zigzag SWCNT could be control 8 . It is shown that 
vacancy could lead to electron spin polarization, hence a 
magnetic SWCNT&iS. But role of doping finite magnetic 
impurity concentration on a zigzag SWCNT is not chal- 
lenged. Our motivation is to investigate possibility of a 
ferromagnetic semiconductor SWCNT in comparison to 
the dilute ferromagnetic semiconductors^. This paper is 
organized as follows: In Sec. II model and treatment for- 
malism of a magnetic impurity doped zigzag SWCNT is 
introduced. In Sec. Ill three cases are considered: in the 
first case just impurity concentration is changed while 
impurity spin magnitude, exchange coupling and tem- 
perature are fixed. We found for the weak exchange cou- 
pling, J, by increasing impurity concentration the semi- 
conductor gap, M, in the spin up density of states is 
closed hence a semiconductor to metal phase transition 
is happened. In the second case just exchange coupling 
is varied and other parameters are fixed, we found by in- 
creasing exchange coupling ,J, the semiconducting gap, 
E' g , of spin up electrons in their density of states is in- 
creased. In the last case just magnitude of spin, S, is 
changed. In this case we found by increasing S the semi- 
conducting gap, EL in the spin up density of states is 
increased. 



Let us consider the Hamiltonian as a general random 
tight-binding model 8 
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where a and (3 refer to the A or B sites inside of the 
graphene Bravais lattice unit cell, where each Bravais 
lattice site is includes two nonequivalent sites that are 
indicated by A and B, c a \ a (cf a ) is the creation (anni- 
hilation) operator of an electron with spin a on Bravais 
lattice site i. and n a — 



iajcr 



c~^ CT c~. is the number operator. 



are the hopping integrals between the 7r orbitals 



of sites i and j with spin a and a respectively, /i is 
the chemical potential, ef is the random on-site energy 
where it takes with probability 1 — c for host sites and 
6 with probability c for impurity sites. Ji is random ex- 
change coupling between 7r electrons and impurity spins 
where it takes with probability 1 — c for host sites and 
J with probability c for impurity sites. Sf is the impu- 
rity spin operator and sf = J2aa' ca \a T ' a& ' c ° ' io' 1S the 
electron spin operator. Where r is Pauli matrix vector 
and Ui is the Coulomb interaction between electrons with 
opposite spin. The matrix form of Eq^is, 

H = - E v " + E *l( v « - rt** + E 

ij i i 

(2) 

where the two-component field operator, \E'] cr , is given by 
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Vj is the random on-site energy matrix operator, 
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and tjj is the hopping matrix in the spinor space defined 
by, 
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and I is a 4 x 4 unitary matrix. 

The electron's equation of motion corresponding to 
Eq|T] in the mean field approximation where the oper- 
ator Si is replaced by it's thermal averaged, Si ~ (Si)thj 
is 

J2 ( (EI - Si + Iih)5u - t a ) G(l,j; E) = I<% (6) 



where G(i,j;E) is the random Green function matrix 
defined by, 
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Here after spin flip is neglected, hence Eq[7| is reduced 
to the two following separated equations for spin up and 
down electrons, 
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(EI - e ia + ml)6u - t iala ) G aa (l, j; £) = Icy 
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where G aa (i,j; E) is a 2 x 2 random Green function ma- 
trix for spin cr defined by, 

c (i a. j?\ — ( G AA (i,j;E) G AB (i,j;E)\ , , 
Gw(y,B) " iG^^i;^ GfBfrjilS))' (9) 

ticrjcr is hoping integral matrix for spin cr which is given 
by, 
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and Ei a is random on site potential matrix for spin cr, 
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where q a = I for spin up and -1 for spin down electrons. 
Eq|Sl could be rewritten in terms of the perfect Green 
function matrix G (i,j;E) as, 

G aa (i,j; E) = G°(i,j; E)+Y / G°(i, I; E)e la G aa (l,j; E) 

(12) 

where G E) is given by 

G°(i,j; E) = ^E e4k ' riJ ( E1 ~ Ck ) _1 ■ ( 13 ) 

k 

with e k = — fil + J2ij t(crjae lk riJ is the band struc- 
ture for perfect system. In our calculations the hopping 
randomness are neglected also we allowed hopping to the 
nearest neighbors and neglected others, so 
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where t ~ 2.5eV is clean system nearest neighbour hop- 
ping integral. Hence 
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and the dispersion relation is 
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where 7(k) = Yli=i e lk ri and are three vectors that 
connect an A (B) site to it's nearest neighbors B(A) sites. 

The Dyson equation for the averaged Green function, 
G aa (i,j;E), corresponding to EalT2lis 

G aa (i,j;E) = G°(i,j; E) 

+ 2 G ( l5 1- E)X aa (l, I; E)G aa (l, 3; E), 
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where the self energy H aa (l, I ; E) is defined by 



(ei a G aa (l,j;E)) = ^'£ l7cr (l,l' ;E)G aa (l' ,j;E). (18) 
l' 

The Fourier transform of G aa (i, j; E) in Ea ll7l is given 
by, 



where 
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S CTCT (k;i?) =-53e- lk ' r «S CT(7 (i,j; J B), (20) 



is the self energy Fourier transform. Since Eas ll9llT7l and 
^Jcould not be solved exactly, we solve these equations in 
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the coherent potential approximation (CPA) . In the CPA 
multiple scattering is neglected and all sites are replaces 
by effective sites except one which is denoted by impurity, 
hence self energy is diagonal Yi aa (i, j; E) = Y l<y<y (E)5ij. 
So Eg 1191 Upl and H5I at impurity site are reduce to, 



0.05 



gu(m;£) = §£(£i 

k 



(21) 



+ G acJ (i,i;E)(e lo 



,{E))GZ P (i,i;E) 
(22) 



and 
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Eas l21l l2*2l and l2*3l are construct a complete set of equa- 
tions that should be solved self consistently to provide 
G„,(i, i; E) and H aa (E). In our calculations the impu- 
rity spin thermal average, (Sf z )th, is evaluated quantum 
mechanically as, 



(sr)t h = 



Eg:: 



, S «z exp (_l J(3 ( n a_ n « )Sr) 



and the average magnetization is calculated from, 
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where £ take all 2 2 impurity configurations and is 
probability of £ configuration. 



III. RESULTS AND DISCUSSION 

To investigate magnetic impurity effects on (10, 0) and 
(20, 0) zigzag SWCNTs, several cases are considered: first 
to obtain impurity concentration effects, average magne- 
tization, M, is calculated in terms of temperature at fixed 
J = — 0.5t, S = t and /i = 0, FigQ]illustrate average mag- 
netization in terms of temperature for different impurity 
concentrations c = 0.005, c = 0.025 and c = 0.05. Also 
in this case at a fixed temperature, T = 0.005i, spin up 
and spin down electrons density of states are compared 
for concentrations c = 0.005, c = 0.025 and c = 0.05. 
Figs[21 and [3] illustrate a semiconductor to metallic phase 
transition in the spin up electrons density of states due to 
increasing impurity concentration for the weak exchange 
coupling J = — 0.5t. In the case of weak exchange cou- 
pling we found by increasing impurity concentration the 
semiconducting energy gap, E^, between valance band 
and conduction band is decreased and a semiconductor 
to metallic phase transition is occurred. While for strong 
coupling the spin up density of states is a semiconductor 
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FIG. 1: Show average magnetization of a (10,0) and a (20,0) 
zigzag SWCNT in terms of temperature at fixed parameters, 
S = 1, J = — 0.5t, /i = 0, u = t and 8 — t, for different 
impurity concentrations, c = 0.005, c = 0.025 and c = 0.05. 



and for spin down is a metal. Second, at a fixed impurity 
concentration, c = 0.025, for variable exchange coupling, 
J = —t, J = —2t and J = —St, average magnetization is 
calculated in terms of temperature. Fig^ shows effect of 
exchange coupling on average magnetization. We found 
for strong impurity couplings, the spin up electrons den- 
sity of states remain a semiconductor with a small gap, 
J57 J , between valance and conduction bands, value of this 
gap is increased for higher exchange couplings , while the 
spin down electrons density of states is continues (it has 
a metallic behavior) . Fig|5] and illustrate our result 
for spin up and down electrons density of states for vari- 
able exchange coupling, J, and fixed other parameter at 
c = 0.025, S = 1, fj, = 0, 5 = t. Third, to observe ef- 
fects of impurity spin magnitude on our system we fixed, 
impurity concentration at, c = 0.025, exchange coupling 
at J = — t, on-site energy at 6 = t and /1 = 0. Magneti- 
zation is calculated in terms of temperature for different 
impurity spin magnitudes 5=1/2, 5=1 and 5 = 3/2. 
FigEI show average magnetization in terms of tempera- 
ture for different spins 5 = 1/2, 5 = 1 and 5 = 3/2. We 
found by increasing impurity spin value, average magne- 
tization is increased but the spin up electrons density of 
states has a gap, E^, which is groves by increasing im- 
purity spin magnitude. Figs|S| and [5] illustrate density of 
states for spin up and spin down electrons in this case. 
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FIG. 2: Show density of states for spin up electrons of a 
(20,0) zigzag SWCNT for different impurity concentrations, 
c = 0.005, c = 0.025 and c = 0.05 at fixed parameters, S = 1, 
J — 0.5t, fj, = 0, w = t and 5 = £. 
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FIG. 3: Show density of states for spin down electrons of a 
(20,0) zigzag SWCNT for different impurity concentrations, 
c = 0.005, c = 0.025 and c = 0.05 at fixed parameters, S = l, 
J = 0.5t, fi — 0, u — t and 5 = £. 



IV. CONCLUSION 

We have investigated possibility of making a ferromag- 
netic semiconductor zigzag SWCNT, by doping magnetic 
impurities such as Fe, Co, Ni, Mn. The doped (20, 0) 
and (10, 0) zigzag SWCNTs are treated in the coherent 
potential approximation (CPA). To find out role of im- 
purity on our systems three cases are considered, first at 
fixed, impurity spin magnitude, 5=1, exchange cou- 
pling at J = — 0.5t and temperature T = 0.005i the im- 
purity concentration are chosen variable. We found by in- 
creasing impurity concentration the semiconducting gap, 
El and El, are closed, hence a semiconductor to metal 
phase transition is take placed. Second, we fixed impu- 
rity concentration at c = 0.025, spin magnitude at S = 1 
and temperature at T = 0.005i, then varying exchange 



coupling, J, we found by increasing J the semiconductor 
gap, E^, in the spin up density of states is increased, 

while gap ,E+, in the spin down density of states dis- 
appeared and have a metallic behavior. Third, at fixed 
impurity concentration c=0.025, J=-t and T=0.005t, the 
magnitude of impurity spin, S, varied, in this case we 
found by increasing spin magnitude , S, the gap in the 
spin up density of state is increased. In summery we 
found by magnetic impurity doping on a zigzag SWCNT, 
one could make a ferromagnetic semiconductor. 
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FIG. 4: Show magnetization of doped (10,0) and (20,0) zigzag 
SWCNTs in terms of temperature at fixed parameters, S = 1, 
c = 0.025, fi = 0, u = t and 8 — t, for different electron 
impurity spin couplings J = — t, — 2t and — 3t. 
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FIG. 5: Shows effects of exchange coupling, J, on spin up 
electron density of states of a (20, 0) zigzag SWCNT at fixed 
parameters, c = 0.025, 5=1, u = 1, 5 = t and /i = 0. By 
increasing J the semiconducting gap is increased. 
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FIG. 6: Shows effects of exchange coupling, J, on spin down 
electron density of states of a (20, 0) zigzag SWCNT at fixed 
parameters, c = 0.025, S = 1, u = 1, 8 = t and fx — 0. By 
increasing J the semiconducting gap ,E\, is completely filled. 
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FIG. 7: Show average magnetization of doped (10,0) and 
(20,0) zigzag SWCNTs in terms of temperature at fixed pa- 
rameters, c = 0.025, J = —t, n = 0, u = t and 5 = t, for 
different impurity spins, S = \ , S = 1 and S = § . 
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FIG. 8: Shows effects of impurity spin magnitude, S, on spin 
up electron density of states of a (20, 0) zigzag SWCNT at 
fixed parameters, c = 0.025, J = t, u = 1, 5 = t and fi = 0. 
By increasing impurity spin magnitude, S, the semiconduct- 
ing gap in the spin up electrons density of states is increased. 
This at high, S, could lead to a phase transition. 
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FIG. 9: Shows effects of impurity spin magnitude, S, on spin 
down electron density of states of a (20, 0) zigzag SWCNT 
where parameters are fixed at, c = 0.025, J = t, u = 1, 
5 — t and /i = 0. By increasing impurity spin magnitude, S, 
the semiconducting gap in the spin down electrons density of 
states is increased, so the system remain semiconductor. 



